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Abstract. We use state-of-the-arts first-principles method to investigate the structural, electronic, and 
magnetic properties of stoichiometric LiFeAs. We optimize fully all the structures, including lattice con- 
stants and internal position parameters, for different magnetic orders. We find the magnetic ground state 
by comparing the total energies among all the possible magnetic orders. Our calculated lattice constants 
and As internal position are in good agreement with experiment. Lhe experimental fact that no magnetic 
phase transition has been observed at finite temperature can be attributed to the tiny inter-layer spin 
coupling. Our results show that stoichiometric LiFeAs has almost the same striped antiferromagnetic spin 
order as other FeAs-based parent compounds and tetragonal FeSe do, which may imply that all Fe-based 
superconductors have the same mechanism of superconductivity. 

PACS. 75.10.-b General theory and models of magnetic ordering, 74.20.-z Theories and models of super- 
conducting state, and 75.30.-m Intrinsic properties of magnetically ordered materials 

1 Introduction a striped antiferromagnetic (SAF) order in the Fe layer 

and a weak antiferromagnetic (AFM) order in the z direc- 
Since the discovery of superconductivity in F-doped LaFeAsC^ 011 perpendicular to the Fe layer. Our calculated internal 
[j], many Fe-based superconductors have been achieved Positions of Li and As are in good agreement with experi- 
by doping appropriate atoms in or applying pressure on ment - Through analyzing electronic and the magnetic re- 
i?OFeAs (R: rare earth) [HI], ^Fe 2 As 2 (A: alkaline earth) sults > we Propose a simple spin model to understand the 

experimental magnetic properties. Thus, we show that sto- 
ichiometric LiFeAs has almost the same SAF spin order as 
other FeAs-based parent compounds and tetragonal FeSe 
do. More detailed results will be presented in the follow- 
ing. 

In next section we shall describe our computational 
details. In the third section we shall optimize fully all the 
structures with different magnetic orders and thus deter- 
mine the magnetic ground state. In the fourth section we 
shall present the electronic structures of the ground state 
phase. Finally, we shall make some discussions and give 
our conclusion in the fifth section. 



2 Computational details 

We investigate all possible magnetic orders for stoichio- 
metric LiFeAs by using full-potential linearized augmented 
plane wave method within the density functional theory 
[55] . as implemented in package WIEN2k 36,37J. The gen- 
eralized gradient approximation is used for the exchange 
and correlation potentials [35]- We treat Li-2s, Fe-3d4s, 
and As-4s4p as valence states, Fe-3p and As-3d as semi- 
core states, and the lower states as the core states. For the 



[SHZHE], SrFeAsF p rTO| f TT V T2] . LiFeAs p ^ll"4]IT5ll8] . and 
■ even FeSe [16111711181115] . The highest transition tempera- 
ture T c has already reached to 55-56 K [4] . Various studies 
have been performed to understand their magnetic orders, 
electronic structures, superconductivity, etc pm f2TIf2"2"ll2"3l 
|24 ] |25 ] |2"6" ] |2"7] . Among the five series of Fe-based supercon- 
ductors, the LiFeAs series seems to be the simplest FeAs- 
based superconductors. Many experiments have been done 
to promote T r [ ^ rT ^ f ^ n^1 J T4 t f28l J 3l l f32 l f33] . It should be 
easier to detect the mechanism of the superconductivity 
in this series. Although there have been some computa- 
tional studies on LiFeAs[8,34.27 , many aspects, even the 
magnetic properties of stoichiometric LiFeAs, are not elu- 
cidated. It is highly desirable to solve such basic issues of 
LiFeAs for further investigations and future applications 
of the FeAs-based materials. 

Here, we investigate the structural, electronic, and mag- 
netic properties of stoichiometric LiFeAs by using state-of- 
the-arts first-principles calculations. We compare all the 
possible magnetic orders, and find the magnetic ground 
state in terms of total energy results and force optimiza- 
tion. Our results show the magnetic ground state features 
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Fig. 1. (color online). Schematic crystal structure of stoichio- 
metric LiFeAs (P4/nmm). The biggest ball is Fe atom, the 
medium As, and the smallest Li. 
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Fig. 2. (color online). The optimization curves of relative to- 
tal energy (meV) as functions of the volume (A 3 ), per for- 
mula unit, for the six possible magnetic orders of stoichiometric 
LiFeAs. 



core states, the full relativistic effect is calculated with ra- 
dial Dirac equations. For valence and semi-core states, the 
relativistic effect is calculated under the scalar approx- 
imation, with spin-orbit interaction being neglected [39]. 
We take R mt x K max =7.5 and make the angular expan- 
sion up to I = 10 in the muffin-tin spheres. We use fOOO 
k points in the first Brillouin zone for different magnetic 
structures. The self-consistent calculations are controlled 
by the charge density in real space. The integration of the 
absolute charge-density difference between two successive 
self-consistent loops, as the convergence standard, is less 
than 0.0001 |e| per unit cell, where e is the electron charge. 
For all the magnetic structures, we optimize the volume, 
geometry, and internal position parameters in terms of 
total energy and the force by the standard 1.5 mRy/a.u. 

3 Fully optimized structures with different 
magnetic orders 

Stoichiometric LiFeAs crystallizes into a tetragonal struc- 
ture with space group P4/nmm at room temperature. As 
is shown In Fig. 1, its unit cell consists of one Fe-As layer 
separated by two Li layers. The experimental lattice con- 
stants arc a = 3. 791 A and c = 6.364A)[14j. To determine 
the ground state of LiFeAs, we arrange the Fe moments 
in the Fe layer so as to form nonmagnetic (NM), ferro- 
magnetic (FM), checkerboard antiferromagnetic (CAF), 
and SAF orders. In the z direction, the successive Fe lay- 
ers can couple ferromagnetically or antiferromagnetically. 
Thus, we have four different AFM orders, namely, CAF- 
FM, CAF- AFM, SAF-FM and SAF- AFM. Structural op- 
timization is done for these six magnetic orders. Our vol- 
ume optimization results are presented in Fig. 2. The equi- 
librium volume is determined by the minimum of the total 
energy against volume. The optimized results summarized 



Table I. The Fe layer in SAF order, either FM or AFM in 
z direction, is lower than that in other magnetic orders. 
We conclude that the magnetic ground state of LiFeAs 
is the magnetic configuration SAF-AFM that the spins 
in the Fe layers are in the SAF order and the inter-layer 
magnetic coupling is a weak AFM interaction. This is the 
same as the magnetic ground states of other FeAs-based 
parent compounds and tetragonal FeSe 8, 12, 18. 22"]l23]. 

The equilibrium volume (88.49A 3 ) of the ground state 
is a little smaller than the experiment value (91.46 A 3 ). 
The magnetic moment per Fe atom is 1.58 /is- The mo- 
ment for the CAF order is smaller. For comparison, we 
do similar calculations using LSDA. The LSDA magnetic 
moment per Fe atom, 0.62/kb and 0.63/is for the SAF-FM 
and SAF-AFM orders, are much smaller than the GGA re- 
sults (1.50/Xb and 1.58/xb), respectively. The Fe-As layer, 
especially the Fe-As distance and the Fe-As-Fe angle, is 
important for the ground state of LiFeAs. This can be 
easily seen in Table 1 that as the magnetic moment M 
of Fe atom increases, the values of the distance d^c-As 
between Fe and As atoms increases and the Fe-As-Fe an- 
gles (a and (i) decreased. The trend of M with different 
dpe-As is consistent with previous result calculated with 
experimental lattice constants [27]. The Fe-As layer as a 
whole is quite robust and Li atoms are dispersed indepen- 
dently between the Fe-As layers, which explains why we 
can get exact positions of As atoms (0.2334 compared to 
the experiment value 0.2365 ), but a little deviated results 
for Li atoms (0.3322 to 0.3459). It looks like that the Li 
atoms play less important roles in the materials. 



4 Spin dependent Electronic structures 

In the following, we study the electronic properties of sto- 
ichiometric LiFeAs in SAF order. The magnetic structure 
in the Fe layer is shown in Fig. 3(a). The arrow implies 
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Table 1. The relative total energy per unit cell (E, with that of the lowest SAF-AFM set to zero), the magnetic moment 
of Fe atom (M), the equilibrium volume per unit cell (V), the internal position parameters of Li and As (z u and z Ab ), the 
distance between Fe and As atoms dp c _As, two angles formed by Fel-As-Fe2 (a) and Fel-As-Fe3 (/3) of LiFeAs in the six possible 
magnetic orders. 



Name 


E (meV) 


M (hb) 




ZLi 


ZAs 


d Fc _As (A) 
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P 


FM 


140. 


0.43 


88.47 


0.3304 


0.2238 


2.345 
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CAF-AFM 


151. 


1.02 
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CAF-FM 
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0.2286 


2.362 


68.3 


105.0 


SAF-FM 


2. 


1.50 


88.31 
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Fig. 3. (a) The magnetic structure of the Fe layer of the 
LiFeAs ground state, (b) The first brillouin zone of the LiFeAs 
ground-state phase with the high-symmetry points labelled. 



that the Fe spins align ferromagnetically in the x direction 
and antiferromagnetically in y direction. Also we show 
the first Brillouin zone with key representative points and 
lines emphasized in Fig. 3(b). The total spin-dependent 
density of states (DOS) and the partial DOSs projected 
in the muffin-tin spheres of Fel, Fe2, Li, and As atoms 
and in the interstitial region are shown in Fig. 4. There 
is no energy gap near the Fermi level and then it shows 
a metal feature. For the SAF order, the spins of Fel and 
Fe2 are antiparallel, it can be easily recognized in DOS. 
The DOS can be divided into two energy ranges: from -6.0 
eV to -3.0 eV and from -3.0 eV to eV (Fermi level). The 
former consists of Fe-3d and As-4p states forming the Fe- 
As bond. In the latter range, the Fe-3d states play a key 
role. It is easily noted that the Li states have very small 
weight in the energy window. The spin exchange splitting 
should be mediated by the d-d direct exchange between 
the nearest Fe atoms and the p-d hybridization between 
Fe and As. 

In Fig. 5 we present the spin-dependent band struc- 
ture of LiFeAs with Fel-3d character in SAF order. The 
plots look like lines, but consist of hollow circles. The cir- 
cle diameter is proportional to the weight of Fel-3d states 
at that k point. Compared with DOS, we can find that 
3d states of Fe play a key role for the magnetic struc- 
ture. It is clear that the bands become narrow along the 
z direction. It means that the electron structure shows 
quite strong two-dimensional character. Also, it can be 
found that there are hole-like (along r~Z) and electron- 
like (along Z-U) band sections. 
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Fig. 4. (color online). Total DOS (states/eV per formula unit) 
and partial DOS projected in the atomic spheres of Fel, Fe2, 
Li, and As and in the interstitial region of the LiFeAs ground- 
state phase. The upper part is for majority spin and the other 
for minority spin. 



The crystal structures with NM, FM, and CAF or- 
ders have the same tetragonal symmetry, while that with 
SAF is orthorhombic. Different symmetries lead to two 
different kinds of crystal field splitting for the Fe-3d state. 
For the tetragonal symmetry, Fe-d state is split into d z 2 , 
d xy , d x 2_ y 2, and d XZiyz , and for orthorhombic symmetry 
it is split into five singlets. From Fig. 5 we can see that 
the bands distributions are symmetric along -T-X and -T- 
Y at -3 eV below the Fermi level, but asymmetric near 
the Fermi level. These are different from the previous NM 
band structures [81 |2~T1[3~4"] ■ in which the band distributions 
are symmetric along _T-X and _T-Y for almost all energy 
range below the Fermi level. 



5 Discussions and conclusion 

The total energy results, as summarized in Fig. 2 and Ta- 
ble I, show that the magnetic interaction in the z direction 
is weak, and The SAF structure is lower by 140-151 meV 
per formula unit than the other magnetic configurations. 
This is consistent with a simple spin model for the Fe 
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Fig. 5. The spin-dependent energy bands of the LiFeAs 
ground-state phase. The left part is for majority-spin and the 
right for minority-spin. The band consists of dots. The bigger 
the dot is, the more Fel-3d character the band at that point 
has. 



spins, H = 52 Uj) JijSi ' Sji where Si is the spin operator 
at site i and Jy coupling constants between sites i and j. 
The summation is limited to the nearest site pairs. is 
J in the y direction, — J in the x direction, and 5 J in the 
z direction. For such spin model, a non-zero 5 is necessary 
to a finite phase transition temperature [401 . The actual S 
must be tiny because no magnetic phase transition is ob- 
served in stoichiometric LiFeAs at any finite temperature. 

In summary, we have investigated the structural, elec- 
tronic, and magnetic properties of stoichiometric LiFeAs 
by using the density-functional-theory method. We deter- 
mine the magnetic ground-state by comparing the total 
energies among all the possible magnetic orders. Our DFT 
results show the magnetic ground state has the SAF order 
in the Fe layer and a weak AFM order in the z direction. 
Our calculated internal positions of Li and As are in good 
agreement with experiment. Through analyzing electronic 
and the magnetic results, we propose a simple spin model 
to understand the experimental magnetic properties. The 
experimental fact that no magnetic phase transition is ob- 
served at finite temperature can be attributed to the tiny 
inter-layer spin coupling. Thus, we show that stoichiomet- 
ric LiFeAs has almost the same SAF spin order as other 
FeAs-based parent compounds and tetragonal FeSe do [5J 
12, 18 , 22 , 23 , which may imply that all Fe-based supercon- 
ductors have the same mechanism of superconductivity. 
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